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Abstract 

We calculate the spectrum of gravitational waves generated during inflation in open 
(Do < 1) inflationary models. In such models an initial epoch of old inflation solves the 
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of the graviton held rapidly approaches the Bunch-Davies vacuum. Then old inhation ends 
by the nucleation of a single bubble, inside of which there is a shortened epoch of slow-roll 
inhation giving Pq < 1 today. In this paper we re-express the Bunch-Davies vacuum for the 
graviton held in terms of the hyperbolic modes inside the bubble and propagate these modes 
forward in time into the present era. We derive the expression for the contribution from 
these gravity waves to the cosmic microwave background anisotropy including the ehect of 
a hnite energy diherence across the bubble wall. [PACS numbers 98.80Bp, 98.80Cq] 
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1. Introduction 


It has been well known for quite some time that gravitational waves are generated during 
inflation and give an observable, and sometimes substantial, contribution to the cosmic mi¬ 
crowave background anisotropy [1-7]. For flat (flo = 1) inflation the spectrum of gravitational 
waves generated and their observational consequences have been studied quite extensively. 
However, for open (flo < 1) inflationary models in which our entire observable universe lies 
inside a single bubble [8,9,10,11], there has been no complete calculation of the gravitational 
waves generated. In this paper we present such a calculation. 

Gravitational waves from inflation result from the stretching of quantum vacuum fluc¬ 
tuations of the linearized graviton held to superhorizon scales. Following a given mode 
of hxed co-moving wavenumber k, one hnds that at early times its physical wavelength 
A = a{t) ■ {27r)/k is much smaller than the Hubble length ijj = H~^ (i.e., the mode is well 
within the horizon). This implies that for determining a physically reasonable vacuum state 
for the mode at early times, one may ignore the expansion of the universe and match onto 
the usual hat Minkowski space vacuum for the graviton held. Once the correct quantum 
vacuum state has been determined at early times, to continue these modes to later times 
become a mathematically well-dehned exercise in classical held theory, which involves prop¬ 
agating the positive frequency modes—those associated with annihilation operators of the 
vacuum—forward in time, through the end of inhation into the present epoch. As a mode 
crosses the horizon its amplitude becomes frozen in. The process of generating gravitational 
waves during inhation is quite analogous to the process of generating scalar huctuations 
during inhation. There is, however, an important diherence. The amplitude of the gravita¬ 
tional waves does not depend on the slope of the potential; rather only the overall height of 
the potential, or equivalently the expansion rate H during inhation, is relevant. As a hrst 
approximation for calculating the gravitational waves, H may be regarded as hxed during 
the relevant epoch of inhation. 

For open inhation identifying the correct initial conditions for the linearized graviton held 
is not as straightforward as for the hat case, because the underlying spacetime geometry is 
more complicated. The simplest case involves one matter held, the inhaton 0, and minimal 
coupling to gravity. More complicated models such as found in [12] have similar properties 
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for the gravity wave calculations done here, so the simple one held model is used in the 
following description. (For a more detailed discussion of single bubble inhation see [10].) 

In the one held models, there is an initial epoch of old inhation during which 0 is stuck in 
a false vacuum with (j) = (j)fy. During this epoch, the spacetime geometry approaches that of 
pure de Sitter space, characterized by an expansion rate Hfy where = {87rG/3)V[(l)fy]. 
During this initial epoch of old inhation the graviton held is driven to the vacuum state. This 
determines the initial conditions using the same considerations as for the hat case described 
above. Then old inhation ends through the nucleation of a bubble, which expands roughly 
at the speed of light. The preferred time slicing inside the light cone of the bubble center 
corresponds to a spatially open universe. Inside the bubble the inhaton held hrst slowly 
rolls down a rather hat part of the potential, giving a shortened epoch of slow-roll inhation 
inside the bubble. Later inside the bubble, the inhaton held rolls more quickly and the 
usual reheating occurs, converting the vacuum energy of the inhaton held into radiation and 
matter. 

The coordinate chart with the line element 


ds^ = —dt^ + a^(t) ■ \d^^ + sinh^[,^] dVL^^2) > 


( 1 , 1 ) 


describes an expanding Friedman Robertson Walker universe with spatially uniform negative 
spatial curvature. Hyperbolic open coordinates (rather than the flat coordinates, to be 
described later) are the natural coordinate choice in the presence of the bubble wall, which 
is why the interior of a bubble is an open universe [8, 9]. For de Sitter space and these 
hyperbolic coordinates, a(t) = sinh[t], and 


ds‘^ 


-dt\ + sinh^[t/i] ■ 


di^ + sinh^)^] 


dVl 


( 2 ) 


( 1 , 2 ) 


where (0 < < +cxd). The hyperbolic coordinate chart has an unphysical coordinate singu¬ 

larity at t = 0, and to determine initial conditions it is necessary to consider a larger region 
of spacetime than that covered by the open coordinates with the line element (1.1). 

The bubble nucleation process underlying open inflation is sketched in Fig. 1, with the 
dashed lines indicating the surfaces on which inhaton held is constant. Roughly speaking, the 


2 



forward light cone of the materialization center M, which we shall call region I and which 
is covered by the coordinate chart just described, may be considered the bubble interior. 
Regions II and III cover the spacetime prior to bubble nucleation and the part of spacetime 
into which the bubble expands, at a speed approaching the speed of light. 



Figure 1—Spacetime Diagram for Open Inflation. Fig. 1(a) shows a spacetime diagram for 
bubble nucleation. The double-dashed vertical line to the left indicates an r = 0 coordinate singularity. 

Time flows upward and the horizontal axis represents a radial coordinate. On the surfaces represented by 
dashed curves the inflaton field is constant. The lower portion of the diagram (with t < 0) represents the 
nucleation of a critical bubble, a classically forbidden Euclidean process. For t > 0 the bubble expands 
classically, at a speed approaching that of light. The classical expanding bubble evolution is SO{3, 1) 
symmetric. In Fig. 1(b) the hyperbolic coordinates that maximally exploit the SO{3, 1) symmetry of the 
expanding bubble solution are sketched. Spacetime is divided into three hyperbolic coordinate patches. 

The light cones separating these regions represent unphysical coordinate singularities of a character similar 
to that of the Schwarzschild horizon. 

For future reference, in Fig. 2, we present the conformal diagram for all of maximally 
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extended de Sitter space. In addition to the regions already mentioned, there also exist 
regions IV and V, which are the past and future light cones of M, the antipodal point of the 
apex of region I. For a discussion of the global structure of the de Sitter vacuum see [13,14], 



Figure 2—Hyperbolic Coordinates for Maximally Extended de Sitter Space. A conformal 
diagram for all of maximally extended de Sitter space is shown here. M is the antipodal point of M. 
The hyperbolic coordinates that exploit the symmetry of the SO{3, 1) subgroup of the full de Sitter group 
SO{4:, 1) that leaves invariant M (and M as well) divides spacetime into the hve indicated coordinate 
patches. 


The subject of scalar perturbations in open inflation has been studied extensively in 
recent years. Lyth and Stewart [15] and Ratra and Peebles [16] calculated the scalar per¬ 
turbations in open inflation assuming conformal boundary conditions for the vacuum of the 
inflaton held as t 0 in region I. Bucher, Goldhaber, and Turok [10] presented the hrst 
computation of the scalar perturbations in open inhation using the Bunch-Davies vacuum 
for the inhaton held outside the bubble during the prior epoch of old inhation as initial con¬ 
ditions and propagating the scalar modes through the bubble wall into the open universe. 
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Yamamoto, Sasaki, and Tanaka obtained essentially the same resnlts nsing Enclidean meth¬ 
ods [17, 18]. Many recent calculations of scalar field perturbations, for example [12,19, 20, 
18, 21 ], take into account additional effects such as finite bubble size, varying bubble wall 
prohle, fluctuations in the bubble wall, et cetera. 

The subject of gravitational waves generated during open inflation was previously in¬ 
vestigated by B. Allen and R. Caldwell in an unpublished manuscript [22], Within their 
approximations they found an infrared divergence in the multipole moments of the cosmic 
microwave background anisotropy (CMB). In their computation, in order to simplify the 
calculation, at early times a flat spacetime geometry (i.e., that of Minkowski space) is as¬ 
sumed, and to improve the infrared behavior of the graviton held during this early epoch, the 
graviton held is given a small mass rUg, which at the end of the calculation taken to approach 
zero. Later on, inside region I on a hyperboloid of constant time with respect to the region I 
coordinates, the spacetime geometry is taken to change discontinuously to that of de Sitter 
space. Subsequently, inhation inside the bubble ends, so that the scale factor 0 ( 77 ) becomes 
governed by a radiation-matter equation of state, and the evolution of the graviton held is 
computed to take into account this change. This is a situation diherent from open inhation 
and violates the Einstein equations for the background solution. 

In section IV of this paper we derive the Bunch-Davies vacuum for the graviton held 
in de Sitter space in terms of the hyperbolic region I coordinates, assuming de Sitter space 
everywhere and a vanishing graviton mass. The Bunch-Davies vacuum for the graviton held 
in de Sitter space, expressed in terms ofthe spatially hat coordinate slicing, is used to express 
the Bunch-Davies vacuum for the graviton held in de Sitter space in terms of the hyperbolic 
coordinates. The latter are the natural coordinates for studying perturbations in an open 
universe. While in principle it should be possible to compute directly the transformation 
between the hat and the hyperbolic modes, in practice this transformation has proven al¬ 
gebraically intractable. Instead, we characterize linear combinations of modes of positive 
frequency with respect to the hat coordinates (which is the same as positive frequency with 
respect to the Bunch-Davies vacuum) in terms of analytic properties in the complex plane. 
In this way linear combinations of hyperbolic modes of purely positive frequency with re¬ 
spect to the Bunch-Davies vacuum may be constructed without explicitly expressing these 
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combinations in terms of the flat modes. The basic approach is somewhat analogous to the 
Euclidean methods used by Sasaki et ah for scalar modes [14], except that here, rather than 
using a Euclidean principle as the starting point, we determine the vacuum in terms of the 
flat modes, using the complex plane as a mathematical tool. 

Our result for the Bunch-Davies vacuum for the graviton field essentially coincides with 
that found by Allen and Caldwell. There are some minor differences: whereas we obtain 
a mixed state directly, Allen and Caldwell obtain a pure state, which in the limit rug 0 
mocks our mixed state. This does not alter the divergence of the CMB multipole moments. 

One caveat in computing the CMB moments relates to gauge fixing. In the Sachs-Wolfe 
formula for the CMB anisotropy the gravity waves contribute only in the integral along the 
line of sight—there is no contribution from the tensor modes from the last scattering surface. 
Since linearized gauge transformations are scalar and vector, one would think that the tensor 
part is gauge invariant.^ However, gauge transformations can move the last scattering surface, 
and thus change the intrinsic contribution in the Sachs-Wolfe formula. If gravity wave 
modes mix with pure gauge modes, cancellations may occur [23]. Linearized gravity wave 
modes typically are taken to satisfy the synchronous gauge condition t^h^p = 0 where t^ 
points along the time direction of the preferred coordinate system. In de Sitter space these 
conditions do not coincide for the flat and for the region I hyperbolic coordinatizations. 
These issues have been investigated for hyperbolic coordinates in flat Minkowski space. 

In this paper we drop the approximation of no energy difference across the wall, which 
is never exactly the case in the presence of a bubble. This requires including hnite critical 
bubble size as well. Taking these effects into account removes the infrared divergence of the 
CMB multipole moments. The calculation resembles the calculations for scalar perturbations 
in refs. [18, 21]. 

The organization of this paper as follows. Sections II and III give the multipole expansion 
for the pure tensor perturbations in hyperbolic space. Section II gives the evolution equation 
for the graviton held, which is solved by separation of variables. In section III we study the 

I We define scalar, vector, and tensor here such that a vector that can be expressed as the spatial gradient 
of a scalar is regarded as a scalar, a tensor that may be expressed as spatial derivatives acting on a 
vector is regarded as a vector, etc. 
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properties of tensor harmonics in three-dimensional hyperbolic space, first written down 
explicitly by Tomita [26]. As mentioned above, in section IV we identify the Bunch-Davies 
vacuum for the graviton held in region I in terms of the hyperbolic modes. In section V we 
compute the effect of nonvanishing bubble size and of nonvanishing energy density difference 
across the bubble wall. This result is used in section VI to give the tensor mode contribution 
to the CMB anisotropy for an open universe, and hnally section VII concludes. There are 
two appendices containing technical details. We set h = SttG = 1 throughout, and H = 1 
until section V. 


2. Gravitational Waves in an Open Universe 

Gravitational waves are huctuations about a background metric. The metric can be 
written as a background metric plus a small perturbation: 

( 2 . 1 ) 

where we set hoo = 0, hoi = 0, = 0, and = 0. These conditions require the scalar and 

vector perturbations to vanish and hx the gauge as well. The unperturbed spatial metric 
(with line element -|- sinh^[^]dG^ 2 )) is used to raise and lower roman indices, and 

the vertical line indicates the covariant derivative induced by 'jij. 

For the background corresponding to the natural coordinates in the interior of the bubble 
center’s light cone, the background metric is given in eqn. (1.2). Hyperbolic conformal time 


V 


In 


tanh[f/i/2] 


( 2 . 2 ) 


will be used primarily in the following, in which the line element is 


ds^ 



—drf -I- d^^ + sinh^ [^] dQ‘^ 2 ) > 


(2.3) 


where —cx) < rj < 0. 
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The condition that the hrst order perturbation of the Ricci tensor vanishes = 0 


gives the equation of motion 


[ 24 ] 




This may be rewritten as 



— 

V(3) 


+ 2/C 


= 0 , 


(2.4) 


(2.5) 


where fC is the spatial curvature, with /C = — 1 for a hyperbolic universe, and we have used 

/ D\ 

^ijkl = ^hik'ljl - Uklil)- Using the ansatz 

7-''(x,r,iC,/’.i. m) = n(C) [T'’4”*(C A C)]/ (2.6) 


we obtain 


vL + (C^ + 3) TU'"(e,^,0;C) = o, 


(3) ' v'’ ' --/J -ij 

<92 + — d ^ + (c^ + 1 ) 

Ct 


ThimO = 0 . 


(2.7) 


The normalization n{C,) will be hxed later. 

Using properties of the hyperbolic tensor harmonics that solve eqn. (2.7)a, the gravita¬ 
tional waves in region I can be expanded as^ 


hrs{v,^,^,<P) = ^ j dC n(C) T/i(r/;C) ap,j,m(C) + h.c. 

jmP 0 

CX) 

= ^ [dC ap,i,m(C) + h.c. 

jmPQ 


( 2 . 8 ) 


Here the operators di{(, P,j,m), aj{(, P,j,m) satisfy the canonical commutation re¬ 
lations [a 7 (C,P,i,m), a/(C',P', = [d\{C, P,j,m), d\{C, P', j',m')] = 0 and 

[hi{C,P,j,m), a|(C',P', = 5(C-C) <^PP' hjf 5mm' where C,C' > 0. 


t It should be pointed out that in addition to the continuous modes with 0 < C < oo there might also 
exist some discrete supercurvature modes, as have been found for the minimally coupled scalar field 
in de Sitter space of mass m when {vP/H'^) < 2. See refs. [14] and [27] for a discussion. The one 
supercurvature mode which would be expected to appear here by direct analogy has zero contribution 
to the CMB because of its lack of time dependence. 



The spatial tensor harmonics and n{() are discussed in the following section. Here the 
time dependence for the hyperbolic tensor modes and the flat tensor modes is found. For 
open de Sitter space, a{ri) = —1/sinh[?7]. Thus for hyperbolic tensor modes eqn. (2.7)b 
becomes 

T" - 2 coth[r/]T' + ((2 + 1)T;, = 0, (2.9) 

replacing —and 1C with + 3) and —1, respectively. Transforming to the dependent 
variable T = sinh^[?7] ■ F, one may recast eqn. (2.9) into the form 

F" + 2 coth[r/]F' + (C^ + 1)F - -^F = 0, (2.10) 

smh [Q 

which is identical to the equation for the spatial hyperbolic radial functions with orbital 
angular momentum I = 1. [See for example ref. [10], eqn. (5.21).] It follows that 

ThivX) = {<sinh[r/] +cosh[r7]} ■ (2.11) 


where ( is allowed to take both signs. 

For the flat tensor modes, which will be needed in order to hnd the Bunch-Davies vacuum, 
the time evolution equation becomes 


rji// rpf 

If - -J- f 


+ — 0 , 


( 2 . 12 ) 


where we replace —and 1C with and 0, respectively. With the substitution T = 

3/2 (+) 

rjj -H, eqn. (2.12) becomes the Bessel equation of order z/ = 3/2 whose solutions Hy^{ujrif) 

and are proportional to the spherical Bessel functions h^^\urif) and \urif) 

1/2 

multiplied by . Consequently, 


Tfinf, w) = /_ . [1 + e '"W 


(2.13) 


where u is allowed to take both signs. 
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3. Hyperbolic Tensor Harmonics 


We now turn to computing and normalizing the hyperbolic tensor harmonics, which 
satisfy the equation 

h.. + (C^ + 3)hy = 0. (3.1) 

The offset in + 3) is chosen for later convenience. The pure tensor character of these 
modes requires that they satisfy the conditions of tracelessness 

K = 0, (3.2) 

and transversality 

= 0. (3.3) 

Here the roman letters {i,j = 1,2,3) indicate spatial indices. 

Since eqns. (3.1)-(3.3) are invariant under rotations and spatial inversion about ^ = 
0, multipole solutions may be classihed according to their angular momentum quantum 
numbers J^, J3, and their parity tt. Parity is either electric with tt = (—)-^ or magnetic 
with TT = (—)'^~''\ denoted hy P = E and P = M, respectively. Fixing j, m, and P, we 
write down the most general symmetric tensor held with these quantum numbers. This 
restricts the angular dependence to a few terms but does not specify the radial dependence. 
Imposing eqns. (3.1)-(3.3) and solving for the .^-dependence gives a solution for each C > 0 
and (j, m, P) for j > 2 unique up to an overall normalization. (There are no monopole 

(j = 0) or dipole (j = 1) modes.) The solution to these conditions has been found by 

Tomita[26]. 

The tensor held with electric parity has the form 

FjiHiC) (ef®e«) 

(3.4) 

+ C) (e“ ® ® e«) VaY,-^(0, y) 

+ T(e;C) (e“®e^) VaV^y,>„(0, y) 

where (a, b = 1,2) indicate angular indices and V indicates an S'^ (rather than an H^) 
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covariant derivative. As a result, (e“ ® e^) = (e^ ® + sin^ 6^ ® e^)- The basis 

functions are = d9, and = d(p. Note that this is not a vielbein (normalized) 

basis. 

We use 


+ 2 coth[,^]Zi)^ + 


sinh^[^] 


+ coi[9]DQ H—. 2 
sin 9 


(3.5) 


and 


= 0 , 

= — coth[,^]e^, 

= — coth[,^]e‘^, 

Doe^ = +sinh[,^] cosh[,^]e^, 
= — coth[,^]e^. 


= — cot[6*]e^, 

D^pe^ = + sinh[,^] cosh[,^] sin^[6']e‘^, 
= + sin[6'] cos6']e‘^, 

= — cot[6']e^ — coth[,^]e^, 


(3.6) 


We now impose the constraints of transversality and tracelessness. For transversality, 
taking the divergence of gives 


V ■ Tjm — 


, o 2coth[^] + 

^ + 2coth|«]F(0 - ^^G(0 - 


j{j + 1) coth[^] 


no 


Yjm{n) X 


sinh^)^] 

^ + 2cothK|i/K) - 


• [j(j +1) -11/(0 


(3.7) 


X 


e dYjm 

09 


+ e' 


ifi dYjm 


d(f 


= 0 . 


Both terms must individually vanish. Likewise, taking the trace and asking it to vanish gives 
the condition 


2G(o jjj + nno 

sinh^[^] sinh^[,^] 


(3.8) 
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Thus transversality and tracelessness give 


m-x) 

mx) 

GjXX) 


sinh^[^] 

j{j + 1) 

sinh^ [,^] 

(J + 2 )(J- 1 ) 
1 






+ 3coth[e]F,-(e;C) 


( dmx) 

V 


+ 2coth[e] 


JU + 0 - sinh () 


mx) 


(3.9) 


The Laplacian in eqn. (3.5) acting on the various components of in eqn. (3.4) gives 


F(OlSm(fi) ■ (e^ ® e^) 

02 p np / n ( n A- 'I'] 

+ 2coth[e]— - ( , 2 J +4coth^[^] ) F(0 


Vsinh^[^] 

[2 cosh^[^]F(,^)Yj'^(r2)] ■ (e^ ® + sin^ 6e^ ® e^) 


YjmX^) ■ {X ® X-) 


(3.10) 


+ [2coth[e]F(0] 


e^ + e^ 


A dYjm , p 

e^)^ + (e' 






d^p 


C^(0^'m(f^) ■ (e^ ® + sin^ ® e*^) 


+ 


4coth^[^] 
_ sinh^)^] 
d‘^G 


GX)Yjm{n) 


(e^ ® e^) 


9^2 

2 cosh[,^] 

sinh3[e] 


2coth[e]^ -2G- 

^ ^ sinh2[e] 


Yjmi^) ■ (e^ ® + sin^ 6e^ ® e*^) 


G(0 


(e^ (g) (g) 

uu 


+ (e^ 




3^) 


dY, 


jm 


dip 


(3.11) 
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H{0 


BY- 

[e^ (g) (g) e^) + (e^ ® e‘^ + ® e^) 

Bo 


BY, 


jm 


Bif 


W±i«y,,„(n)//K).(e^ 0 e^) 

smh [Q 

B‘^H 


B^ 


„ - 2H -Acoth^inH - t 
2 sinh 2 [e] . 


X 


(e^ ® ® e^)-^r' + ® e^ + e‘^ ® 

Bo Bip 


+ 4i7(0 coth[e] 


^ 0^d‘^Yjm , , f 9‘^Yjm , . ^ ^dYjm 


\ (96*(9v? 


(e^ ® + e^ ® e^) (- cotO^ 


B(p 


(3.12) 


and 


B(p 


V2 /(OI (e^ ® e^) + (e^ ® + sin 9 cos 9^^ 

+ (e^ ® e‘^ + e^ ® e^) 

\ B9Bip 

^ -2/(0,(, + l)coth^a (^) (e,^e^) 

smh [^] 

+ ■ (e^ ® e^ + sin2 9 e‘^ (^ e^) 

smh [^] 

^ 2[j(j + 1) - 1] coth[^]/(0 
sinh2[^] 


X (e^ ® e® + e^ ® e^) + (e^ 


B9 


BY- 

(e^ ® ® e^)—^ 

Bip 


X 


L/l/ L/^y 

^ - 2 coth[e]^ + - 2/(0 coth 2[0 - 

[de sinh2[0 sinh2[0 

^ fB%m , . ^ ^BYjm 

)^ + (e^ ® eO + sm ^ cos 


YTY^O 


(e^ Y) 


(e^ ® e‘^ + e^ ® e^) (- cot 9^^^ 


(3.13) 


To solve for T) (0 0) we take the (e^®eO component of the Laplacian acting on eqn. (3.4), 
and after applying the substitutions in eqn. (3.9), the coefficient of ® term in (3.1) be- 
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comes 


£!!5^ + 6coth|5]^?S^+ [(C^ + 3) + 6coth2|5]-^^^^]F,«;C) = 0. (3.14) 

c's c's L smh [Q _ 

With the change of variable 4 >j{^]C) = sinh^[,^] ^), one recovers the differential equation 

for the sca/ar hyperbolic radial functions (see, for example, [10] eqn. (5.21)): 

^1^00 + ^ ^ 1 ) _ Ml _ 0 , ,3,15) 


Consequently, 


sinh [^] 


• n ( 1 -^ 


(3.16) 


where the normalization Nj{() is determined in the following, and we have imposed regularity 
at the origin. 

For future reference, for j = 2 


P(c.n- 

and from eqn. (3.9) 


3C^coth[^]cos[C^] + (C^ + C-3Ccoth2[^])sin[C^] . (3.17) 


G'2(^;C) = ^2(0-sinh)^] ■ _2^2 

+{ "f' ~ (c^ - 2 + 3 coth^Kl) }cos|Cfl . 


^^ 2 (^;C) = iV 2 (C)• 


^2(^;C) = ^ 2(0 ■ sinh[^] ■ 


^■[{V(C^ + 4-6coth2[e])}cos[Ce] 

{^^^^^^(2coth2[^]-(2 _2j|sin[C^] , 
[fC^cothie]/ ^ 


(3.18) 


-5 - 2C^ + 3 coth^)^]) |cos[C^] 


-C^-4C3-3C , Ccoth^ie] 


(^- coth2[^] + (2 ^ 2^ I sin[C^] 
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Similarly, for the magnetic parity, the tensor held must take the form 


9, ip; C) = () (e“ ® + e? ® e“) LaYjm{e, if) 

(3.19) 

+ C) (e“ ® ® e“) 

The magnetic parity modes do not contribute to the CMB anisotropy because their compo¬ 
nent along the (e^ ® e^) direction vanishes; therefore, we do not give their explicit form. 

Normalization of the Hyperbolic Tensor Harmonics. To normalize the tensor har¬ 
monics, we impose the condition 


oo 

J J dn^l^)Tij{x;C,P,j,m)* T*^(x; C', P', j', m') 

0 

= d(( — (') ■ Sp^p/ ■ 6jj' ■ 6m,m' 


(3.20) 


where P indicates mode parity and x = (,^, 6 , 0). 

Because the tensor harmonics are eigenfunctions of a self-adjoint operator, the inner 
product is proportional to a ^-function. Thus Nj(() is determined by the coefficients of 
in the asymptotic expansion of the tensor harmonics for large (For 0, Fj, Hj, Gj, Ij — > 
0.) In comparing the asymptotic behaviors of F, G, H and I, it is more meaningful to 
consider the rescaled quantities Fj(C;C) = ^j(^;C), = G'j(C; C)/smh^)^], = 

0 /sinh[,^], and sinh^[,^], components with respect to a normalized 

‘vielbein’ basis. From eqn. (3.9), F sinh ^[^], H ~ sinh ^[^], and G, I sinh ^[^]. 
For ^ ^ X, G and I dominate, exactly as one would expect. That is, at large distances a 
spherical gravitational wave should locally resemble a plane gravitational wave propagating 
in the radial direction. 

To compute the coefficient of the (5-function in eqn. (3.20), we impose a boundary con¬ 
dition at .^ = ^max (for specihcity say Dirichlet boundary conditions) and take the limit 
^max OO- For ^max ^ 1; fhe integral in eqn. (3.20) is dominated by the G and / com¬ 
ponents. These may be approximated by their large-.^ asymptotic forms, starting with 
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eqn. (3.16) and substituting sinh[,^] —(e^/2) = {w/2). This gives 


F,(e;C)^4iV,-(C) 



[yj+iC + yj-iC^ 


= 4 Nj{() e ■ {iC)j + c.c. 


(3.21) 


Here (x)j is shorthand for x(x — 1)... (x — j). Using eqn. (3.9), we obtain 


C) ^ ^ ■ [«) ■ ] > 


j(j + l) 

2 Nj(0 e+^ 


^7(C;C) ~ ^-r ■ [« + 1) ■ «) • (O,- + c.c. ' , 

^ j(j + l)(j-l)(j + 2) J 

0 ^ ■ [« + 1) ■ «) ■ (^0J + C.C. ■ . 


(3.22) 


Consequently, for (.^ 3> 1), 


T(x; (, E,j, m) ^ / -r ■ (^C + 1) ' (^C) ' {K)j ' 


X V~aVlYjm{e,<t>) + + l)Yjm{e,4>) ( e“ ® 6^ 


+ C.C. 


(3.23) 


Inserting this asymptotic expression into eqn. (3.20) and factoring out 5jj' 5mm' gives 


sinh^)^] / dfi Ty(x;C,U, j,m)T*^’(x;C,^, j,m) 




2 4C"(C" + 1) [C^(C^ + U)...(C^ + J^)] 

j2(j + l)2(j _ l)2('j + 2)2 


(3.24) 


dQ 57jYjm + T 
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For the angular integration, note that 


dQiVaV^^Yjm) = - / dQiV~^Va'^r.Y, 


a'/ b^jm 


dn[v~aV~^^Yj^ 


(3.25) 


= -J dn(Vb{V^ +-RjYj^) 

= [ j(j +1) ] ■ b'(j +1) -1 ], 


given that [Vg, V^] Vg/ = g°‘'^R~b/Vdf = Ri{V^f) = \R V^/. Since the unit two- 
sphere is isotropic, = \d^^R. Also, for the two-sphere R = 2. Thus the integral on the 
last line of eqn. (3.24) is equal to [j{j + l)][^j{j -l- 1) — 1], and eqn. (3.24) reduces to 

\AT.(/-\\‘2 _ /-2(/-2 I -[\\/-2(/-2 I t2'i (/-2 | ^■2'ji ^max 2g^ 


from which it follows that 

^ _ Vd{j + l)(j ~ l)(j + 2) 

^ CVi + VC^(C^ + i2)... (C^ + j2) ■ 


(3.27) 


The TT comes from passing to the continuum and the phase of Nj{() has been chosen to 
be real. 

Flat Tensor Harmonics. The general form of the flat modes is needed to identify the 
Bunch-Davies vacuum in the next section. These may be regarded as the large-^, small-.^ 
limit of the hyperbolic modes, a limit in which the effects of spatial curvature disappear. We 
use lower case letters to denote the flat analogues of hyperbolic quantities. In particular. 


f{r,u;) = nj{u;)-Y U- — 


cos cur 


(3.28) 


Similarly, eqn. (3.9) is modihed to 




(j + 2)(i-l) [“ <9r ' r 

9 j{r;uj) = - [j{j + l)ij{r-) - r^fj ]. 


^ dhj 4 

2 ^r- + -hj - fj , 


(3.29) 
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In particular, for j = 2 


f2{r;uj) 


g2{r-uj) 


/i2(r;a;) 


nir-^u) 






cos [cur] + 


cos [cur] + 


-3cn 

+ 32 ) 




Q Sen 1 . r I 

cn H-^ sin [cur J, 


-o;2 o;4 


2r2 6 


H-cos knr + 


-u 


2r 


cos[cc;r] + 


-^5^2 

12 


-UJ'^ OJ 

_ + _|smKrl. 

u 

- TT + TT^ I smlcnri. 

6 2 r^ 


(3.30) 


In order to calculate the normalization, one needs the asymptotic behavior for large r: 


^ nj{C) ■ ^ + c.c. 


hj{i]u) ^ n,(C) ■ -TT 


^nj(C) ■ -T- 


i(i + l) 2 rL 
1 


(io;)2+2 6+*^"^ + c.c. 




j(j + l)(j-l)(j + 2) 
1 


■ r- 


2 (j-l)(j + 2) 


■ r- 


(icc)2+3 6+*^"^ + C.C. 


•^)j+3 g+it.r 


(3.31) 


Again ij and gj dominate, and following the same steps as for the hyperbolic modes we get 


nj{uj) 


Vj(j + i)(j + 2)(j-i) 


(3.32) 


This is the C —^ 3> j limit of eqn. (3.27). 

Normalization of the Time Dependent Mode Functions. We may dehne the anti¬ 
symmetric bilinear form 


(W, V) = - y W«/3 (E) (^^p)V“^(S), (3.33) 

s 

where 7/, V are solutions to eqn. (2.4). This product is analogous to the Klein-Gordon 
product for scalar held modes. Eqn. (2.4) insures that the current ^/^^(S) (iVp)V"^^(S) 
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is conserved and thus that (W, V) is invariant under deformations of the surface S. In order 
that the modes for a spacetime M orthonormalized with respect to (3.33) are associated 
with operators that satisfy the customary canonical commutation relations, it is necessary 
to choose S in (3.33) so that S is a Cauchy surface for M. The product ( , ) is the same 
for all Cauchy surfaces for the spacetime M. A Cauchy surface is a spacelike hypersurface 
which each non-spacelike curve intersects once and only once [28]. In this paper we shall 
consider both the case where M is just region I, in which case a surface of constant region I 
hyperbolic time serves as a convenient choice of Cauchy surface, and also the case where M 
is all of maximally extended de Sitter space, in which case the surface defined by r = 0 in 
the region II hyperbolic coordinates serves as a convenient Cauchy surface. 

For the hyperbolic modes dehned in eqn. (2.6) in region I one calculates 


(T(C,P,j,m),T(C',P',/,m')*) 

= - f d^x a^{r]) 7;^(x,r/;C,P, j,m)(iSQ) 7j(x,r/;C',P', j',m') 


(3.34) 


= ''^^(C) 2 C (C^ + 1) ■ <^(C - O ■ ^P,P' ■ ^j,j' ■ ^rn,m' 


where Pq = (l/a)P^ = — sinh[?7]P^. Note that eqn. (3.34) is independent of rj. This may be 
seen by applying Drj to eqn. (3.34) and using eqn. (2.9). For the mixed tensor representation 


chosen above, the covariant derivative Pq may be replaced with the ordinary derivative d^. 
So, reading off. 


«(C) 


1 

V2C(C2 + 1)‘ 


(3.35) 


As the initial conditions are determined by the bubble which extends outside of region 
I, a proper Cauchy surface for initial conditions for an open universe extends outside of 
region I as well. It was shown in [14] that for some cases, inner products taken on a Cauchy 
surface in region II agreed with those on hxed time surfaces in region I and V, even though 
the latter do not make up a Cauchy surface for the whole spacetime. The norms agreed 
for scalar helds with sufficiently fast falloff at inhnity, a condition satished by subcurvature 
modes of the Laplacian (modes with eigenvalue + 1 > 1). There are some differences 
in extending this comparison between norms taken in I and V and norms taken in II for 
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gravity waves. When the gravity waves are continued across the light cone, into region II, 
time and space are interchanged and the gauge choice becomes hufi = 0 instead of hrj^ = 0. 
As M is a spacelike coordinate, this is not the usual gauge for metric perturbations. It may 
be thought of of as analogous to axial gauge (where A3 = 0 ) rather than Coulomb gauge 
(where Aq = 0) in electrodynamics. Secondly, the inner product in region II involves a 
wronskian in (the analytic continuation of) and is more complicated due to the tensor 
structure in As shown in appendix B, this inner product on a Cauchy surface in region II 
coincides with the inner product in regions I, V, for modes which have sufficiently fast falloff 
at inhnity. The helds with sufficiently fast falloff for gravity waves are again subcurvature 
modes, 0 < < cx). As the falloff for the gravity wave modes is as fast as that for the 

scalars in regions I, V (both go as sinh~^[^] for large .^), it was reasonable to expect this. 


4. Identifying the Gravitational Wave Bnnch-Davies Vacuum 


The preceding section gave the mode expansion for the linearized gravitational waves 
in region I hyperbolic coordinates, which are the natural coordinates for the open universe 
inside the bubble. In this section the initial Bnnch-Davies vacuum is expressed in terms 
of these open hyperbolic modes. In the open universe inflationary scenario, the Bnnch- 
Davies vacuum is a preferred quantum state for de Sitter space in the sense that it is a 
weak attractor: any initial quantum state for perturbations from de Sitter space, subject 
only to the requirement that the initial energy density be hnite, approaches the Bnnch- 
Davies vacuum to arbitrary accuracy after a sufficient amount of inflationary expansion. 
The convergence is weak rather than strong because the initial perturbations are not erased 
but rather pushed to larger and larger scales, so that for an observer able to probe only a 
hxed physical volume, the perturbations seem to disappear. 


The Bnnch-Davies vacuum is physically characterized using the flat coordinates for de 
Sitter space. In these coordinates, at sufficiently early times, a mode of hxed co-moving 
wavenumber evolves as if it were a mode in Minkowski space. For the hat coordinates the 
line element is 





+ n 


2t 


f . 


dr'j + rj dQ^ 2 ) 


(4.1) 
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where —cx) < tf < +cx), or in terms of flat conformal time rjf = —e 


-if 


ds^ 




—drj'j + dr‘j + rj 


(4.2) 


Early on, when the mode is subhorizon, so that there are many oscillations within an 
expansion time, one identihes the mode with positive frequency asymptotic behavior with 
an annihilation operator of the Bunch-Davies vacuum. These are the modes that behave as 
^-tkrif ioY Tjf —CX), where rjf is flat conformal time. 



Figure 3—Flat Coordinates for de Sitter Space. The region covered by the flat coordinates is 
shown in a conformal diagram for all of maximally extended de Sitter space, identical to that in Fig. 2. 
Although the flat coordinates cover only half of maximally extended de Sitter space, the diagonal line, 
which is the boundary of the region covered by the flat coordinates, represents an initial value surface for 
all of de Sitter space. 


To relate the flat and hyperbolic coordinates one can embed de Sitter space in (4+1)- 
dimensional Minkowski space as described in appendix A. Using sinh[t/j] = —1/sinh[?7]. 
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cosh[t/i] = — coth[? 7 ], and coth[t/j] = cosh[ 17 ], one gets that the relation between the coordinate 
systems is 




rif 


sinh[^] 

cosh[,^] + cosh [ 77 ] ’ 
sinh [ 77 ] 

cosh [ 77 ] + cosh[^] 


(4.3) 


The flat coordinates do not cover all of de Sitter space but only half of maximally extended 
de Sitter space, as indicated in the conformal diagram in Fig. 3. Nevertheless the flat 
coordinates cover enough of de Sitter space to contain an initial value surface for all of de 
Sitter space. The null surface indicated by the dashed diagonal line, the boundary of the 
region covered by the flat coordinates, is such a surface.^ 


To identify the hyperbolic modes of positive frequency with respect to the Bunch-Davies 
vacuum, it is convenient to use the null coordinates in the flat chart 


Uf = 7]f + r/, 
^f=Vf-rf. 


(4.4) 


One has nj < 0 and \uf \ < —Vf. We also define the region I hyperbolic null coordinates 


U = V + ^: 
v = v-^ 


(4.5) 


where < 0 and \U\ < —V. Using the above relations between and ( 777 -,rj) one can 

see that two sets of null coordinates are related according to 


U 


In 


1 +^/ 

1-Uf 


V 


In 


'1 + W ' 
-l-W- ' 


(4.6) 


Conversely, Uf = tanh^ and Vf = tanh ^ . Thus we see that the region I hyperbolic 
coordinates only cover the region —1 < Vf < 0, \uf \ < —vj in terms of the flat coordinates. 


f Technically, because this surface is null it is not a Cauchy surface, but initial data on this surface can 
be continued everywhere in de Sitter space. 
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The positive frequency modes with respect to open hyperbolic light cone coordinates in 
region I are 


e 


iCU 


1 +^/ 
1-Uf 


-iC. 


1 


-KV 


' 1 + ^/ ' 


(4.7) 


However, a full mode function for the Bunch-Davies vacuum must be specified in both regions 
I and V. A region I mode function must be continued into region V in order to calculate its 
inner product on the full space and the choice of analytic continuation distinguishes between 
positive and negative frequency Bunch-Davies modes. The identification of the positive 
frequency modes for the Bunch-Davies vacuum is determined by the analytic properties of 
the factors in and appearing in the hyperbolic mode functions above. The other 

factors in the hyperbolic mode functions contain only isolated singularities and no branch 
cuts, and thus will be seen to be irrelevant in making this identideation. The bubble interior 
lies within the strip —1 < Uf < -fl. Outside this strip one encounters a branch cut, taken 
here to lie near the real axis starting at Uf = 1, passing through -|-cx), coming back from 
—oo, and finally ending at Uf = —1, as indicated in Fig. 4. Thus [(1 + Uf)/{1 — w/)]*^ has 
two possible continuations to the entire real line—one through the upper half-plane above 
the branch cut, and another through the lower half-plane below the branch cut. 
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Figure 4—Analytic Continuation of the Mode Functions. The points uf = — 1 and uj — +1 

are branch points. To impose single-valuedness, we take a branch cut to start uj = -|-1, extend to 
= -1-00, come back from uf = -1-00, and finally end at uf = — 1. The analytic continuation through 
the upper half-plane onto the remainder of the real line corresponds to the Bunch-Davies positive frequency 
modes. Similarly, the continuation through the lower half-plane corresponds to the Bunch-Davies negative 
frequency modes. 


One can identify the positive frequency Bunch Davies modes by requiring they have 
vanishing overlap with the flat negative frequency modes (cr > 0). Thus one requires 

= («) 

— CxD 

for all td > 0. This is satisfied by the analytic continuation through the upper half-plane 
above the branch cut, as deforming the contour toward +ioo makes this integral vanish. 
As a result, positive frequency Bunch Davies modes correspond to analytic continuation 
through the upper half plane into region V. This identihcation generalizes to gravitational 
wave hyperbolic modes and other more complicated mode functions of the form 

(l^) (l^) 

where R{uf,Vf) and S{uf,Vf) are rational functions. 


-I-CO 

!■ 




Figure 5—Dual Pair of Flat Coordinates. Both panels are conformal diagrams for maximally 
extended de Sitter space. The shaded triangle in Fig. 5(a) shows the region covered by the flat null 
coordinates (ttj, vf), while in Fig. 5(b) the region covered by null coordinates (uj, vf), is shown. 
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We shall use two sets of flat coordinates {uf,Vf) and {uf,Vf), as indicated in Fig. 5, 
related by a spatial reflection symmetry of de Sitter space that maps region I into region V. 
It may be shown that 

Uf = —, Vf = —, (4.10) 

W «/ 

either directly by using closed coordinates or by the analytic methods in appendix A. 

This transformation relates the flat coordinates {uf, Vf) covering the upper left triangular 
wedge of de Sitter space to the coordinates {uf, Vf) covering the upper right triangular wedge, 
as indicated in Fig. 5. As discussed earlier, for positive frequency mode functions the analytic 
continuation is with positive imaginary part, (—1) = e*^, and so 


l-UfJ 

l-VfJ 




( 4 . 11 ) 


As a result, positive hyperbolic frequency modes in region I become negative 

hyperbolic frequency modes in region V with an attentuation factor and negative 

frequency region I hyperbolic modes become positive frequency hyperbolic modes in region 
V with an amplihcation factor The result of all this is that the positive frequency 
modes in region I are the modes proportional to for all values of ( plus the 

continuation into region V using the upper half plane. The negative frequency modes are 
the same collection of modes with the opposite analytic continuation into region V. 

We may normalize the modes by combining the fixed time surfaces of regions I and V. 
As shown in [14] for scalar subcurvature modes, and in the appendix for gravity waves, for 
the subcurvature modes under consideration here, the inner product is equivalent to that 
taken on a Cauchy surface for all of de Sitter space. Let P, j,m) denote the modes 

with positive hyperbolic frequency in region I, as dehned in eqn. (2.6), normalized using only 
the region I hxed time surface, and let be the analogously dehned modes in 


25 



region V. It follows that in region I one has the mode expansion, 


Pjm 0 


dc 


+ 1) 

e-C /2 




3 + '^C/2 


s[(e+^^ - e-<)V2 


ThivX) bi{C,P,j,m) 


{e+< - e-^C)V 2 
and similarly in region V 


Th{v,-C) bv{C,P,j,m) 


+ h.c.} 


(4.12) 


^ f d( 

t:l vW+i) 




3+'^C/2 


3+7rC 


g-7rC)l/2 


ThimC) bv{C,P,j,m) 


e-C /2 

(e+< - e-<)i/2 


Th{v,-C) bi{C,P,j,m) 


+ h.c.} 


(4.13) 

where the annihilation operators bi{(, P, j, m) and bv{C, P, j, m) annihilate the Bunch-Davies 
vacuum for the graviton field and satisfy the usual commutation relations 


[bA{C,P,j,m), bB{C',P',j',m')] =0, 

= 0, (4.14) 

|i>/l(Ci4’> j>m)> i\{C.',P','j\ rn')] = Sab 4(C “ CO Spp' Spy Smm' 


where {A, B = I, V). In the Bunch-Davies vacuum we observe a doubling in the number of 
modes. This reflects the existence of correlations between regions I and V—or, alternatively, 
between the inside and the outside of the bubble. 

Another way of understanding this result for the vacuum is to note that the scalar mode 
functions and mixed tensor (with one raised and one lowered index, T^) mode functions 
considered here have the same time dependence. As the choice of vacuum depends only 
on the time dependence of the mode functions, the time dependence in the Bunch-Davies 
vacuum also has the same form for both cases.* 

For the scalars there is a discrete vacuum mode [14] and here one also is normalizable, 
by the same arguments. 

-k J.D.C. thanks K. Schleich for discussions about this. 
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5. Gravity Waves and a Bubble with Finite Energy Difference 


In this section we take into account the effects of the nonvanishing size of the critical 
bubble and of the nonvanishing difference in energy density across the bubble wall. In the 
presence of a bubble, one starts with a Bunch-Davies vacuum with H = Hp as an initial 
condition outside the bubble and then continues these modes across the bubble wall into 
the open universe, where we idealize the expansion rate to be a constant with H = Hp. Of 
course, in realistic single-bubble inflationary models H inside the bubble is neither constant, 
nor does it correspond to the true vacuum as the subscript T suggests. However, in order to 
capture the qualitative consequences of changing H across the bubble wall without getting 
involved in the messy details of specific models, we consider the idealization of an infinitely 
thin wall, with constant H = Hp inside, for calculating the gravitational waves. Whereas 
earlier in the paper we set H to unity, from here on factors of H will be displayed explicitly. 
Given that the gravity waves and the scalars have the same time dependence, and this is 
the only place the dependence on the wall appears, one can carry over the results for scalar 
fields. In the case at hand, we have a massless held before and after tunneling and a finite 
energy difference across the wall. The specific relevant results are from [18, 29, 21] and are 
briefly sketched below. 


As will be shown below, any bubble with Hp ^ Hp has a finite radius and extends outside 
the light cone (region I) covered by the (? 7 , 6, 0) coordinate system. By considering a curve 

of constant held value outside the light cone in Fig. la, one sees that the matching across 
the wall thus has to be done in region II and its Euclidean continuation. The continuation 
into region II from region I is —>■ —^ r + inj2, and so the metric in region II 
becomes 


ds‘^ 



du^ — dr^ + cosh^ 


(5.1) 


For an open universe a{u) = l/(hf cosh)^]). The analytic continuation of the basis functions 
is given in somewhat more detail in appendix B and is straightforward. With the Euclidean 
continuation appropriate for tunneling, one also, before nucleation, rotates the space to 
Euclidean time r —> irp. 
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The wall has the same description in the Euclidean and Lorentzian regions in spacetime, 
since it only depends on u and not r, te- For a thin wall, the u dependent scale factor is [30] 


a{u) = 


Ht cosh[u +51 ’ for « > UR, 


(5.2) 


where 


Hr COsh[M + 5] 


He 


cosh[M — Ur] cosh UR + sinh[M — ur] ^cosh^[Mi^] — {Hr/HeY 


The metrics are explicitly matched at the same value of a{u), the bubble radius 

R= ^ ^ 


He cosh[-Ui^] Hr cosh[Mj^ + 5] 


(5.3) 


(5.4) 


The mode functions in this background have only their -u-dependence changing across 
the wall. Requiring continuity in the unchanged (r, 9, 0) dependence forces the exterior 
mode functions to match onto interior modes with the same value of To match the 
mode functions across this wall, we start with the inital u dependence corresponding to the 
analytic continuation of the time ( 77 ) dependence found in eqn. ( 2 . 11 ). As the exterior of the 
bubble has H = Hr, the false vacuum mode functions ^), extend to region II as 


A^(m) = He cosh[M] ^tanh[M] — iC, 


(5.5) 


up to a factor of i. These are matched across the wall onto modes of the true vacuum (as 
there is some slow roll after tunneling this is not exactly the true vacuum but is close enough 
for our purposes). The corresponding mode functions interior to the bubble, with u > ur, 
H = Hr, are 

Bi^{u + (5) = Hr cosh[M + 5] ^tanh[M + 5] — • (5-6) 

Matching these and their first derivatives at the wall, at a{uR) = R, gives 


A(^{u) = ai^B(^{u + 5 ) + (3(^B_(^{u + S) 


(5.7) 
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where 


(5.8) 

13 ^ = J_e*C(2««+<5) 

Zl(^ 


and 

2 ; = tanh[M7^] — tanh[M7^ + 5] = y^l — {Hj?R)‘^ — y^l — {HtRY- (5.9) 


Note that these become trivial in the limit (5—^0, that is, when the energy difference 
{Hp — Hp) —^ 0. From the dehnition of 2 :, it appears possible that 2 : = 0 is possible even 
ii Hp ^ Hp, via taking R = 0. However, within the new thin wall approximation of [10], 
R oc with a nonzero constant of proportionality, and in the usual thin wall 

approximation [31] one has that the radius obeys 


R^ = 


s! 


Co(/f2 _ H^) + ci^2(^2 + ^2) + 


(5.10) 


where Si is the surface tension (related to the integral of the change in the background held 
through the wall) and cq, ci, C2 are positive constants. Since both Hp and Hp are hnite and 
the surface tension does not vanish, R is hnite. 

One actually has to go to another coordinate system in order to show that there is a ‘time’ 
where there is only false vacuum and no bubble [30] for initial conditions. The consequence 
here is that the mode functions must be orthonormalized once they are continued across the 
bubble wall, as they do not correspond to a normalized initial mode functions in this other 
‘time.’ Thus the mode functions given in eqn. (5.7) are not yet properly orthonormalized. At 
the nucleation time (conventionally taken as the r = 0 slice in region II), the whole system is 
rotated to Lorentzian signature. The mode functions can then be orthonormalized in region 
II, or since these are subcurvature modes, one can continue them to regions I and V and 
orthonormalize them there. The latter was done in detail in [18, 29] and the former was 
done in detail in [21]. Substituting the tensor rather than scalar spatial functions into these 
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results, in region II the wavefunction is, (again up to overall factors of i), 


jm i 

X<{ bi{C,EJ,m) ■ 


dC 


V4<XC^Tiysinh< L 

+ *S'(+C) \ B(^{u) 


+ i7r/2,0,0,r/;C) 


■{c(+C);3< + S(+C) a_<;} 




-0 a-c + ^(-C)/9c 5-c(«) 


+ C(-C) /9-c + ^(-C) «c 5c(«) 


(5.11) 


As a result, the positive frequency part of the wavefunction inside region I for subcurvature 
modes is 




dC 


jm 0 


V2C(C2 + 1) 




X< bi{C,E,j,m) ■ 


■) g+'>TC/2 

“c + >^(+0 / 5 - c } _ ^- n (: y /2 +0 

/3c + ^(+0 «-c} ^^+^C-e-<)V2 


+bv{C,E,j,m) 


-0 «_c + ^(-C) 


,-</2 


+ |C(-C)ff-c + 5(-C)ac| (,.,c_,-,C)l/2 


(e+^C - e-^C)i/2 

3 + '^C/2 


Tft(r/; -0 


Tfe(r/;+C) 


(5.12) 


where 

Ei{C) = - [kcl^ + = Di{-Q = 1 + 

D2{C,) = a^(3^ = ^{2tC,-z)^<^-, 


(5.13) 
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and 


c^(C) 




-^(0 = 


-D2 


Di 


m \i di-\d2\^ \ 


1 - 4 1 - 


\M 

Dl 


-D2 

m ^ 


1 - 


4C' 


;^2 + 4^2 


(5.14) 

The identity D\ — |/l 2 p = Di was used in the above and can be checked by substituting 
into the dehnitions. 


For scalars, an analogue of the vacuum supercurvature mode in the presence of a bubble, 
with (C = —i, remains normalizable when matched across the wall!^*^ This is true for the 
case here as well, but as this mode is rj independent, it does not contribute to the CMB 
anisotropy below. Consequently we do not quote its value here. 


6. Implications for the CMB Anisotropy 

Gravitational waves provide a time dependent background for the photons in the CMB 
[32] (for a review see [33]). The contribution of gravitational waves to the CMB anisotropy 
is given by the integral in the Sachs-Wolfe formula [32]: 

Vo 

= ^J *1 = 1)0 - n, n) 

\ a 

” 5 /*' 

Ve 

The metric h is dehned with the conformal factor scaled out, hence h^i, = h^^/a^irj) in our 
notation. Here we have chosen the path to be parameterized by conformal time, where rjQ 
is the observing time and rjf, is the last scattering time for the photon. The surface term 
vanishes for the tensor contribution and hence is omitted. The radial-radial component of 
the mode function has spatial dependence proportional to as seen in eqn. (3.4). 
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The Bunch-Davies positive frequency part of the temperature contrast operator is, using 
eqn. (5.12) from the last section, 
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2 J 

jm I 


1 

+ 1 ) 



ho -h;C) 


X 


bl{C,E,j,m) 
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Thim -0 
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( 6 . 2 ) 


where the dots indicate derivative with respect to conformal time. The Bunch-Davies nega¬ 
tive frequency component is 


5T^gw\^) 

Tcmb 


Tcmb 


(6.3) 


It is customary to expand the CMB anisotropy in terms of multipoles according to 


STg^n) 
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(6.4) 


The statistical average of the ensemble of classical gravity waves is found by taking the 
corresponding quantum average, so that the two-point correlation is 
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(6.5) 

for all /. In showing this, it is useful to note that |C(C)P + |>S'(C)P = 1 and the dehnitions 
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in eqn. (5.13). Also since Di = ^[1 + |q;^P + |/9(P] and — |/3(^p = 1, 


1 + l«c^ + l/^cP 


< + 1 . 


( 6 . 6 ) 


Consequently as one probes larger wavenumbers C, (through larger-^ multipoles) the influence 
of the bubble dynamics quickly becomes negligible in a uniform manner. In other words, 
effects of the wall radius and finite energy difference rapidly disappear below the curvature 
scale, and thus are expected to be hidden in the cosmic variance, just as in the case for the 
scalar fluctuations. 


The integrand in the eqn. (6.5) for the CMB multipole moments is well behaved for small 
C,. This can be demonstrated with the explicit form of T2('Cj C) (3.17). [Recall that for 

gravity waves the 7 = 0, 1 moments vanish because the graviton is a spin-two particle.] One 
then notes that Fj for higher j are obtained by taking derivatives, which will not alter the 
leading power of C (although the coefficients may change). As C —^ 0, the time dependence 
factor in the basis functions becomes 


dr/Tfii^hC)] — 9^[{K sinh[j/] + cosh[j/])e i/rsinhl?/] (1 - + 0(C^ 


( 6 . 7 ) 


The first term in the curly brackets has the form 
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since 77e[l — i((r]i — 772 )] = 1. For the second term in the curly brackets, for small (, 
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The argument in brackets on the right has the form 
Tie 
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including the small-^ expansion of the exponentials. Thus, after factoring out the sinh^[? 7 ] 
dependence, the term in curly brackets in eqn. (6.5) behaves as 


1 ^ + 0 ( 1 ) 
irC(l + |4) 


( 6 . 11 ) 


for Q small. In the Sachs-Wolfe integral above, we also have the factor of in the measure. 
Since Nj ~ one has (using eqn. (3.17) in conjunction with eqn. (3.27) ) Fj ~ constant. 
This implies that the integrand approaches a constant as C 0, rendering the integral 
infrared convergent. 


As {Hp — Hp) becomes small, both and 2 : approach zero for hxed C, and the sec¬ 
ond term in curly brackets approaches zero. For vanishing ; 2 , as is found in the vacuum, 
eqn. (6.11) approaches coth[7rC], making the integrand appear to diverge as ~ We do 
not have an intuitive understanding of this limiting behavior. Vanishing is never the case 
in the presence of the bubble since it implies exactly zero energy difference across the bubble 
wall. 


In order to calculate the airm oae needs the time dependence of the wavefunctions from 
the inflationary period, through radiation domination and into the current phase of matter 
domination. Matching conditions have been found by [22] and are calculated analogously to 
the flat [1-7] and closed [34] universe cases. Time dependence has also been considered in 
[35] but seemingly for a different initial vacuum. 


7. Discussion 

We have determined the initial condition for the graviton held in an open universe 
originating from a bubble inhation model and calculated the contribution from gravitational 
waves to the CMB anisotropy. The total observed CMB anisotropy for a given multipole 
is obtained by combining in quadrature the gravity wave contribution calculated here with 
the scalar held contributions for the particular model. The ehects of the bubble wall for the 
tensors, just as for the scalars, seems to be conhned mostly to very large scales, corresponding 
to C small. It appears that the ehects of curvature provide a much larger ehect for the tensors 
than the ehects of the wall!^**' 
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Noted Added: After this work was completed we learned that Allen and Caldwell and 
Sasaki et ah have reached similar conclusions with respect to the influence of nonvanishing 
bubble size. 
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APPENDIX A—Relation Between Flat and Hyperbolic Coordinates 

For the flat coordinates the line element is (setting H = 1) 


ds^ = —dt^j -I- ■ dr'j -|- r'j dC?^2) 


(A.l) 


where —oo < tf < -f-cxD, or in terms of flat conformal time rjf = —e R 



(A.2) 


To relate the flat and hyperbolic coordinates one can embed de Sitter space in (4-1-1)- 
dimensional Minkowski space (see ref. [10], section 5 or ref. [28]). The Minkowski coordinates 
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are {w, u, x, y, z) = (ta, m, r) and de Sitter space is defined by 


- 2,-2 -2 1 
r + M — w =1. 


(A.3) 


The embedding of the open hyperbolic coordinates is 

w = sinh[t/i] cosh[,^], 

u = cosh[t/i], (A.4) 

f = sinh[t/i] sinh[,^]. 

Generally, f > 0, and as 0 < < cx) and 0 < .^ < cx), so that one sees that region I 

hyperbolic coordinates cover the range 0<tc<cx, 1 < u. The fiat coordinates are 
embedded in (4+l)-dimensional Minkowski space according to 

-1 f 


tf = ln[t(; + m], Vf = 


— -V, Tf = —-T- 

w + u w + u 


(A.5) 


These cover w,u > 0, a larger region than the hyperbolic open coordinates. As sinh[t/i] = 
— 1/ sinh[7f], cosh[t/i] = — coth[? 7 ], and coth[t/j] = cosh[7f], the relation between the coordinate 
systems is 

sinh[^]_ 

fnl ’ 

(A.6) 

Sinn 17/1 

Vf = 


T f = - 

■' cosh [^] + cosh [if] ’ 
sinh [if] 


cosh [if] + cosh[,^] 

The fiat coordinates do not cover all of de Sitter space but only half of maximally extended 
de Sitter space, as indicated in the conformal diagram in Fig. 3 in the text. 

The null coordinates in the fiat chart are 


(A.7) 


Uf=7]f + r;, 

^f=Vf-rf. 

One has nj < 0 and \uf\ < \vf\. We also define the region I hyperbolic null coordinates 

U = V + ^, 

V = y-i 

where V < 0 and |17| < \V\. Using the above relations between (if,0 {j}f,rf) one can 


(A.8) 
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show that the two sets of null coordinates are related according to 


U 


In 


1 + ^/ 
l-Uf 


V 


In 


1 + ^^/ 
1-Vf 


(A.9) 


Conversely, Uf = tanh ^ and Vf = tanh ^ • Thus we see that the hyperbolic coordinates 
cover only the region \uf \ , |n/|<l, \uf\<\vf\m terms of the flat coordinates. 

In order to see the continuation into region V, consider the transformation u —u in 
the embedding Minkowski space. This can be accomplished by taking th = —ty — In 
terms of conformal time in region V, this is 


Tjv = In tanh(—t/j/2 + m/2) = ln(—1) + In tanh(t/j/2 — m/2) 

= In(-l) - ijj 


(A.IO) 


Thus we have 

U = Tjv + ry = In(-l) - V 
V = rjy — ry = ln(—1) — U 


(A.ll) 


where U, V are the hyperbolic open coordinates in region V. The dehnition of ln(—1) requires 
a choice of analytic continuation, which has been identified in terms of Uf,Vf, so converting 
to these coordinates, one has 


Uf = tanh(f//2) = tanh(—C/2 — ■iim/2) = — coth(t//2) = —Vj ^ 
Vf = tanh(C/2) = tanh(—C/2 — ±m/2) = — coth(f7/2) = —uj^ 


(A.12) 


as given in the text. This transformation relates the flat coordinates {uf,Vf) covering the 
upper left triangular wedge of de Sitter space to the coordinates {uf,Vf) covering the upper 
right triangular wedge, as indicated in Fig. 5. 
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APPENDIX B—Inner Product in Region II 


For this appendix, H = 1. To continue into region II, > —ie > r + i7r/2, and 

the metric becomes 


ds^ = a^{u) du^ — dr^ + cosh^[r] dP 


( 2 ) 


(B.l) 


where a{u) = 1/ cosh[M]. The Wronskian, with respect to —iVr on the symmetric perturba¬ 
tions of the metric, gives for an inner product 


dudncosh^lr] ^ ^ ^ ^ ^ ■/ c 

-- 2 —-W^(r, S;C,J,m) [idr) (r, S;C ,J ,m) 

cosh \u\ 


(B.2) 


The Bunch-Davies vacuum mode expansion [ see eqn. (4.12) ] in region I 
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pjml y2c(c^+i) '■ 
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{e+< - e-<)V2 


Th{v,-0 bv{C,P,j,m) 


+ h.c. 


becomes in region II 


r r p-|-7rC/2 /'A'l 
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PJm g ^ ^ 2 
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(•BA) 

Note that the complex conjugate is only taken in the u dependence in region II and both 
terms have the same r dependence corresponding to positive frequency. The integral over 
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the surface S in region II is an integral over the analytic continuation of time u and over 
(6*, 0). The analytic continuation of Tii{r]) multiplies the whole expression and is independent 
of r. So for calculating the Wronskian, first consider only the r dependence and the integral 
over 6, 0. The integral over u will be done subsequently. We have 

= Fj{T + iTT/ 2 ;C) (e^ 0 e^) YjmiO,(p) 

+ Gj{T + i7r/2; () 5~^i (e“ ® e^) 

+ Hj{T + in/2; () (e“ ® ® e“) VaYjmiO, 99 ) 

+ Ij{T + in/2; C) (e“ ® e^) <^) 

Thus 

T; = -F,-(r + *7r/2;C)lS'm 
TJ =-i/,(r+ *7r/2;C) 

T“ = cosh-2[r] ^ C) V“Fjr„ 

T0 = cosh-2[r] (J,-(r + *7r/2; C) + *7r/2; C) 


(B.5) 


(B.6) 


As we have V‘^Yjrn = —j{j + l)h^m and from earlier 

j dn{yaViYjmy = -J dfi(v“VaV^F,>n)* (vV,w) 

= -y' dfl(VaV^VW,>n)* (vV,w) 

= - J dfl(v^|v2 + iF}Fj>„)* (v%m') 

= [ i(i + !)]•[ i(i + 1 ) - 1 , 


we obtain immediately that 


j dVtT^rT^* = {Fj^rF* - 2j{j + 1) cosh-\]Hj*^rH* 

+ cosh ^[T](/j d rljUij + l))(j(j + 1) ~ 1) + 2Gj d rGj 
-j{j + l){Gj*^rIj + 

= {iFj'SrP; - 2j(j + 1) coslr^lTl//,-^,//; 

+ i cosh""*- 2)ljl}rIj }^jfSmrn'- 


(B.7) 


(B,8) 


The sign in front of Hj is due to the negative signature of r. There is a factor of two in front 
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of both Hj and Gj. For Hj there is both a term Tia as well as a term T^i, and for Gj there 
is a contribution from both G 2 and G 3 . The minus sign in front of the last term is due to 
the negative sign relating and 

After some algebra, one hnds that 


dnr^'d'rT^ 




2C"(l+C")cosh-‘|T] 
j(j + l)(j + 2)(j-l) 


FjdrF’SjfK 


(B^9) 


After some work one can show that 


FjdrF; = Nj{C) 


2 *sinh[ 7 rC] .3,, , , ^2^ ,.2 , ^2^ 


cosh°[r 


-c'>(i+c")( 4 +c")'"(i"+n. 


(B40) 


Putting this all together and substituting for we get 


rrB* 2iCsinh[7rC] 


dn T§ d rTr = 


TT cosh^ [r 


2 r 


(B.ll) 


The integral over u remains, with mode functions 


n(C)ef 


, . -u^ ^ urn iCu cosh[n] n(C) 

Th{\n{-te );C) / . . = «-tanh[n]) e*- ==, 


(B.12) 


so the last integral, including the rest of the measure, becomes 

j' ^ cosh^[r] cosh[n] (< - tanh[M])n(C) cosh[M] (-<' - tanh[M])n(C') 
J cosh^ [u] a/2 sinh[7rC] a/2 sinh/C'] 

27r cosh^ [r] 

4 Csinh[7 rC]‘ 


(B.13) 


As a result, we have 


- j dticosh ^lu]msh^\T]U^{u,E;C,j,m) (Af r)U^‘(u, ,j' ,m') = - (')SjfSmm' 

(B44) 


which agrees with eqn. (3.34). 
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